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Abstract 

It has recently been shown that the neutrino can have non-zero electric 



> 

fT^ ■ charge in a number of gauge theories, including the Minimal Standard Model. 

Cp ■ Assuming non-zero neutrino charge, we develop a new approach to the solar 



neutrino problem. The key idea is that the charged neutrinos will be deflected 
by the Lorentz force while they are crossing the solar magnetic fields. Such a 



Oh! 

^ • deflection will result in the anisotropy of the solar neutrino flux. Because of 



this anisotropy, the solar neutrino flux registered on earth can be reduced as 
compared to the Standard Solar Model prediction. The mechanism is purely 
classical and does not require neutrino oscillations, spin-flip or neutrino decay. 
We discuss qualitatively the consequences of our scenario for present and fu- 
ture solar neutrino experiments as well as differences between our mechanism 
and other proposed solutions. 
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I. INTRODUCTION 

A long-standing puzzle at the crossroads of elementary particle physics and astrophysics 
is a deficiency of the fiux of solar neutrinos registered in four different experiments. The data 
of the two-decade-long CI — Ar experiment at Homestake |jl[ have more recently been com- 
plemented by the direct observation of solar neutrinos at the Kamiokande water Cherenkov 
detector ||^ as well as by two different Ga — Ge experiments, SAGE ||^ and GALLEX [Q. 
Being sensitive to different parts of the solar neutrino spectrum, these four experiments 
together provide the data which are difficult to understand with a single set of solar param- 
eters chosen within the Standard Solar Models |^,||. Chances are that we have to modify 
some aspects of neutrino behavior in order to explain the observations. 

An incomplete list of candidate solutions proposed so far includes neutrino oscillations 
possibly enhanced by Mikheev-Smirnov- Wolfenstein mechanism [^§, neutrino decay [Q, 
neutrino spin-precession ||10[ in the solar magnetic field and resonant spin-fiavour conversion 



scenario 



jrT|. Recently, a more sophisticated version of the last scenario has been pro- 
duced which employs more detailed assumptions about the small scale structure of the solar 
magnetic field in the convective zone |1^ . 



In this work we would like to analyze a simpler and perhaps more natural hypothesis 
TB| in connection with the solar neutrino problem: the possible existence of nonzero electric 



charge of the neutrino. In fact, the neutrino is the only elementary particle, besides the gauge 
bosons, whose electric charge is normally assumed to be zero. But if the neutrality of the 
gauge bosons is deeply rooted in the principle of gauge invariance, there are no compelling 
reasons whatsoever for the neutrino to have zero charge. 

Of course, the neutrino is assumed to be exactly neutral within the Standard Model. 
However, the recently^ developed approach to the problem of the electric charge quantization 



^Early attempts to study the possibility and consequences of charge dequantization date back 



to Einstein |14| who noted that a tiny difference between the electron and proton charges could 



has led to the reahzation of the fact that in a fairly large class of gauge models, including the 
Minimal Standard Model, the electric charge can be dequantized [|T^ (see also ^^). This 



means that the electric charges of elementary particles can take different values from those 
conventionally assumed: Qu = 0, Qi = — e, Qu,c,b = 2e/3 and Qd,s,b = — e/3 (e being the 
modulus of the electronic charge). In particular, the neutrino can acquire nonzero electric 
charge. (Another interesting aspect of the theories with dequantized electric charges is that 
one might speculate about the possibility of time dependence of the electric charges within 
such theories |1T8[| .) 

To be more specific, consider the Standard Model with the electric charge dequantized 
through the formula [0 

Q = Q standard + ^{Le — L^j. (1) 

Within this model, the electron and muon neutrinos acquire electric charges of the same 
absolute value but of opposite signs, while the tau neutrino remains neutral: 

g(z/e) = -Qiu^) = ee; g(z/,) = 0. (2) 

Another possibility is to dequantize charge by a similar formula fl^ 



Q — Q standard + c(-t>e ~ -t>rj- ("jj 

In this case, z/g and z/,- are (oppositely) charged while z/^ is neutral. (We would like to 
emphasize that the above examples are by no means exhaust all possible ways of giving the 
neutrino nonzero charge.) 



In other words, in Refs. |16| it has been shown that the Standard Model contains an 



additional free parameter e which must be determined experimentally along with the other 



account for the Earth's and Sun's magnetic fields. Later on, Lyttleton and Bondi ||T^ proposed 
the idea that such a tiny difference could cause the expansion of the Universe due to electrostatic 
repulsion. It took experimentalists several decades to rule out these hypotheses (and some of their 
modifications). 



more familiar parameters such as Higgs mass or Yukawa couplings. Of course if it were found 
that e is nonzero, it should be very small anyway (see next Section) and that would create 
one more hierarchy problem. Yet taking into account the existence of a few such problems 
already, the appearance of a new one does not seem strong enough argument to disregard 
the possibility of nonzero e. 

Furthermore, one might argue that nonzero neutrino charge does not follow from any 
theoretical principle, whether established or hypothetical (with the exception of the well- 

^ ^ we 

know that the zero neutrino charge does not follow from anywhere, too! 

Another possible objection against particles with small fractional charge is that it is 



known rule "all which is not forbidden is allowed"). But now, based on the works [16 



difficult to embed them into grand unified theories |T9|. Yet theories with paraphoton 
provide a viable alternative p^ . 



So, at present the cases of zero/nonzero neutrino charges must be considered as two 
working hypotheses on the equal footing, only experiment being capable to provide the 
ultimate answer. The situation with neutrino charge is very similar to the situation with 
neutrino mass: while zero mass is the prediction of the minimal standard model, most 
physicists agree that the question of zero/nonzero neutrino mass has much more to do 
with experimenting than with model- building. While it would not be easy to detect the 
neutrino charge, the consequences of such discovery should certainly be dramatic, ranging 
from the prospects of detecting relic neutrino through its electromagnetic interaction to 
possible better ways of managing neutrino beams, creation of neutrino optics etc. 

Finally, let us note that in the present work we are not concerned if the neutrino mass is 
zero or not. Certainly, there exist well-known difficulties associated with charged massless 
particles |21| . However, one can take a pragmatic point of view [^ and keep developing a 
theory until one runs into any inconsistency. No such inconsistency seems to show up in our 
treatment. An alternative point of view is to give the neutrino a Dirac mass by introducing 
additional Higgs multiplets. 

Note also that even if the neutrino is massless in vacuum, it cannot be considered massless 



inside plasma. This is because the vacuum dispersion relation i? = |p| gets changed by the 
weak interaction of neutrino with plasma |p . In other words, there arises a refraction index 
for the neutrino propagating through plasma. Thus, the situation with infrared divergencies 
might be better for a neutrino in plasma than in vacuum. 

To conclude, it seems that assuming nonzero neutrino charge is certainly not more heresy 
than assuming nonzero neutrino magnetic moment, or mass and mixing angles. 

In addition, there exist quite an independent motivation to study the behaviour of a 
charged neutrino inside the Sun. The point is that the neutrino electromagnetic properties 
get modified by plasma effects and under certain conditions these modifications result in 



an induced electric charge of the neutrino [^ . We stress that it happens in the Minimal 



Standard Model where the neutrino has zero intrinsic electric charge. 

The purpose of this paper is to draw attention to the fact that the possible existence of 
a very small electric charge of the neutrino may be a clue to the solar neutrino problem^. 
The idea is that the charged neutrinos are deflected by the solar magnetic field while passing 
inside the Sun. Thus the resulting neutrino flux is made anisotropic which leads to the solar 
neutrino deficit registered on the Earth. (Within our scenario, this deficit is not real but 
only apparent in the sense that the total Att solar neutrino flux is not changed as compared 
with the standard solar models.) 

Section 2 summarizes the existing bounds on the electron neutrino charge. In Section 
3 we give an order-of- magnitude estimate of the effect. Section 4 is devoted to the more 
detailed calculation of the solar neutrino deficit in the present scenario. In Section 5 some 
experimental implications are discussed and one of them - the prediction of a second neutrino 
flux - is briefly outlined in Section 6. Our main results are summarized in Section 7. 



^The possible role of the charged neutrino interaction with the terrestrial electric field, in con- 
nection with the solar neutrino problem, was discussed previously by G.C.Joshi and R.R.Volkas 



(unpublished). 



II. PHENOMENOLOGICAL LIMITS ON THE ELECTRIC CHARGE OF THE 

ELECTRON NEUTRINO 



For a recent detailed analysis of various bounds on the neutrino charges see 
The strongest model- independent^ constraints on the electron neutrino charge, Q(z/e 
ee, come from three sources. 



1) Analysis of the data on v^e elastic scattering ||22|| : 

e < 3 X 10"^°. (4) 

2) Study of the electromagnetic neutrino production in the core of the Sun and its 
influence on the solar energetics [^ : 

e < lO^^l (5) 

This estimate follows from the fact that if the neutrino is charged then the neutrino- 
antineutrino pairs will be produced in the decays of plasmons (which can be viewed as 
massive photons) in the core of the Sun. The neutrino will then escape freely from the Sun 
thus taking away a lot of energy which would be inadmissible in terms of the solar energy 
balance. 

3) Analysis of the data on SN1987A supernova explosion P3[: 



These upper limits on the electron neutrino charge e were based on the experimental 



detection [^] of the neutrino signal from the supernova 1987A explosion. It was claimed 
that if e were larger than (^ then the intergalactic and galactic magnetic fields would 
lenghthen the neutrino paths, and neutrinos of different energy could not arrive on the 
Earth within a few seconds of each other, even if emitted simultaneously by the supernova. 



^By model- independent we mean the constraints that do not rely on additional assumptions such 
as charge conservation or the equality Q{y^ = Qiy^). 
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Yet it is generally believed (see, e.g. [^) that the constraint based on SN1987A ar- 
guments, although stronger, is less reliable than the previous ones because it involves the 
details of the galactic magnetic field which are not very well known. 

There exist even more severe, but less direct constraints. They are based on the experi- 
mental data on the neutrality of atoms and neutrality of the neutron. These data give limits 
on the sum of the proton and electron charges [^: Q{p) + Q{e) = (0.8 ± 0.8) x lO^^^e and 
the neutron charge [^: Q{n) = (—0.4 ± 1.1) x 10~^^e. Then, assuming charge conserva- 
tion in the neutron beta-decay n — > p + e~ + z7g we can obtain the bound on the electron 
animeutrino charge: Q^C'e) < 3 x lO^^^e. Finally, assuming validity of CPT symmetry^ with 
respect to z/g and z7g charges, one can claim that 

g(z/e) < 3 X lO^^^e. (7) 

In particular, this constraint is valid within the model with charge dequantization through 
Eq. (|l|) and Eq. (^). Yet the requirements of the electric charge conservation and CPT 
symmetry, although very general and perfectly valid up to now, are themselves a subject of 
experimental testing at present^. Furthermore, there exist several models in which charged 
neutrino arises as a natural consequence of the electric charge violation WM (in those specific 



^Besides charge conservation and CPT, a number of usually unspoken but very important assump- 
tions underlying the constraint (0) are made. For instance, one has to assume that the electric 
charges of free electrons and protons are exactly the same as those of atom-hound electrons and 
protons. Another fundamental assumption, as noted in Ref. |^], is that the electric charge, as 
measured by interaction with an electromagnetic field, coincides with the electric charge assigned 
by the charge-conservation law (see also a discussion of that point in [^). According to Ref. p^ , 
it is possible to construct models in which it is not the case. Under ordinary circumstances there 
is no doubt in the correctness of the above axioms, but when it comes to such fantastic accuracies 
as 10^^^, it does not seem unreasonable to question those axioms, too. 

^Note that it is possible to constrain the Ve charge assuming only charge conservation in the decay 



models the constraint (|^) is still not invalidated). 

In this work we are not concerned with the problem of constructing a viable model which 
would allow one to avoid the constraint (|^). Instead, in the main part of the paper we do not 
use any model- dependent arguments and we treat the neutrino charge e as a free parameter; 
also we consider some interesting consequences of the choice e ~ 10^^^. Of course, should 
any of the above constraints change in the future, our results would be easily reformulated 
by simply rescaling the magnitude of e. 

A massive Dirac z/g with an electric charge of 10~^'^e and a mass of < 7^9 eV [^ would 
have a Dirac magnetic moment yU > 5 x 10^^ fiB, where fiB is the Bohr magneton. Therefore 
one might suspect a contradiction with the existing upper bounds on the neutrino magnetic 



moment obtained recently by a number of authors (see, e.g., [^ and references therein): 
//(i^e) < 10^'-^'^^^^'' Hb- However, these limits apply only to the anomalous magnetic moment, 
but not to the Dirac magnetic moment. In our case, the neutrino anomalous magnetic 
moment due to the electromagnetic radiative correction is equal to e^ajlv: and therefore 
negligible. Still, some of the constraints on /i(z/e) can be translated into limits on Qiv^^. 
For example, the analyses of plasmon decay into neutrino-antineutrino pair due to /i(t'e) in 
various astrophysical contexts can be carried out also for the case of a charged neutrino. 
Yet, the resulting limits are not significantly different from the solar constraint (|]). 

Note also that recently there has been considerable interest in discussing the possible 
existence of new particles carrying very small electric charge (" milli-charged particles") 
||r?| . These works contain detailed discussion of many phenomenological constraints on such 
particles obtained from a variety of sources (including astrophysics, cosmology, geophysics 
and macroscopic electrodynamics). Many of those constraints apply to the case of electron 
neutrino, too; we shall not repeat that material here. 



n — 5- p + e"*" + z^e, but not the equality Q(y^ = Q{ve). Naturally, this constraint turns out to be 
much weaker than (^), namely: Q{v'e) < 4 x lO^^e [^. 



III. SIMPLE ESTIMATE 

Since the neutrino charge, if any, must be so tiny, ^ lO^^^e, at first sight it seems unhkely 
that it would play any role at all. To make a rough check whether it is true or not, let us 
start by a dimensional order-of- magnitude estimate of the effect of the solar magnetic fieldQ 
on the the charge neutrino motion inside the Sun. As a dimensionless characteristic of the 
effect it is natural to choose the neutrino deflection angle. The deflection angle 7 may be 
written as the product of the deflection rate, dj/dt, and the time of flight through the 
convective zone, r: 

7 ~ — - X r, r ~ — . (8) 

' dt c ^ ^ 

where D ~ 2 x lO^^cm is the thickness of the convective zone, c is the speed of light. Now, 
in a crude approximation, the value of the deflection rate d'j/dt will depend only on the 
neutrino electric charge, ee, the neutrino energy i? ~ 1 MeV and some characteristic value 
of the solar magnetic field. He- The only combination of dimension time~^ which can be 
made out of these quantities is the Larmour frequency, 

eecHc , , 

uJl = — ^, (9) 

(We use the Gauss system of units throughout the paper, i.e., e = 4.8 x 10^^° esu, 1 MeV 
= 1.6 X 10~^ erg.) Putting all together and adopting e ~ 10"^^ for the neutrino charge and 
He — 10^ -7- lO^G for the solar magnetic fieldQ we finally obtain for the deflection angle 

7 ^~ 0.01 ^0.1. (10) 

E 



^Here, by "the solar magnetic field" we mean the large scale toroidal magnetic field in the con- 
vective zone. For the moment, we ignore the magnetic fields in the core and the radiative zone - 
not because they are irrelevant but because we do not know much about them. 
This choice is discussed in more detail below, between Eqs. (30) and (|3l]). 



Thus we see, that contrary to a naive expectation, the possible tiny charge of the neutrino 
may indeed play some role in the neutrino propagation through the solar magnetic field. 
Yet from our estimate it is hard to see the relation between the neutrino defiection angle 
and the neutrino deficit observed on the Earth. Also, the above consideration suggests that 
the seasonal variations of the solar neutrino fiux must be large (because of two facts: first, 
the Earth's orbit is inclined with respect to the solar equator; second, the solar magnetic 
field is small at the equator and it grows away from the equator). It will turn out that this 
suggestion is wrong and that seasonal variations in our scenario would not be significant. 
So it is worthwhile to study the problem more carefully. 

IV. CALCULATION OF THE NEUTRINO FLUX ON THE EARTH 

We have seen, by an order-of-magnitude estimate, that the suggested value of the neu- 
trino charge , e ~ 10^^'^, may be relevant to the solar neutrino problem. Now we can start 
a more exact calculation of the reduction of the solar neutrino fiux at Earth due to the 
bending of neutrino trajectories by the solar magnetic field. The key quantity we need to 
find is the angular distribution of the neutrino fiux coming out from the surface of the Sun 
assuming that the initial flux coming out from the central region of the Sun is isotropic. 

Therefore, the procedure is three step. First, we start with the isotropic angular distri- 
bution. Second, we find out how the neutrino trajectories are bent by the solar magnetic 
field. Mathematically, this bending of trajectories can be viewed as a change of variables 
in the distribution function. The last step, then, is to make explicitely such change of an- 
gular variables that would lead us to the final angular distribution. From this anisotropic 
distribution we shall easily obtain our key result, Eq. p^) : the reduction of the neutrino 
fiux predicted by the present hypothesis (and then compare it with the reduction of the fiux 
observed experimentally). 

For definiteness and simplicity of the calculation, we make the following assumptions. 
First, we assume that all the neutrinos are emitted from a point source located in the centre 



10 



of the Sun. In other words we neglect the radius of the neutrino-emitting zone, R^,, as 
compared to the solar radius Rq (their ratio is R^/Rq ~ 0.1). Next, we suppose that the 
neutrino deflection angle is small (the precise meaning of that assumption will be discussed 
below). Furthermore, what we are interested in is not the complete 47r angular distribution 
of the neutrino flux but only a small part of it within the angular interval swept by the line of 
vision Earth-Sun. This is a narrow interval close to the solar equatorial plane: —7° < 6 < 7°. 
This fact will be used throughout the calculation to simplify it. 

Finally, a few remarks about the solar magnetic field. The structure of the solar magnetic 
field is very complicated and is probably one of the worst known aspects of the Sun. Yet, 
some important information about it can be obtained from the observations of solar activity 
such as sunspots which strongly depends on the behaviour of the large- scale magnetic field 
inside the Sun. For instance, observing Zeeman effect for the light emitted from the sunspots 
one can measure the strength of the magnetic field near the surface of the Sun which can 
then be translated into estimates of the asimuthal magnetic field in the convective zone of 
the Sun. In what follows we will rely only on relatively well understood features of the solar 
magnetic fields which include (see, e.g. ^^): 

l)the existence of the large-scale toroidal magnetic field in the convective zone of the Sun 
(toroidal means the magnetic lines of force are closed circles parallel to the Sun's equatorial 
plane) ; 

2)the directions of the toroidal fields in the nothern and southern hemispheres of the Sun 
are opposite. Besides, these directions reverse themselves every 11 years; 

3)the strength of the toroidal magnetic field is nearly zero in the equatorial plane and 
it grows with the separation from the equatorial plane within the narrow transition zone, 
between about +10° and —10° latitudes. Therefore, in this tranzition zone there exists a 
gradient of the toroidal magnetic field in the direction parallel to the rotation axis of the Sun. 
Unfortunately, the magnitude of the gradient is not very well known and we will consider it 
as a parameter varying in a certain range (see below, between Eqs. ( pOD and (PT|). It is this 
gradient rather than magnetic field itself that plays the key role in reducing the neutrino 
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flux on the Earth within the present scenario. 

According to our plan, we start with the isotropic neutrino flux from the centre of the 
Sun described by the angular distribution 

— = No = const, (11) 

or equivalently, 

where 6 and are the usual angles in the spherical coordinate frame whose origin is placed 
in the centre of the Sun, y-axis points to the Earth and z-axis is the rotation axis of the 
Sun. After passing the convective zone fllled with magnetic fleld the neutrino is deflected 
from the direction characterized by the pair of spherical angles 6 and to the new direction 
characterized by the spherical angles a and (3. Using the equation of motion, we can ex- 
plicitely write down the angles a and (3 as functions of 9 and 0. Conversely, we can express 
a and (3 as functions of 9 and 0: 

9 = T{a,f3),(f) = P{a,(3). (13) 

Now, all we need to get the new distribution function from the old one is to change variables 
from 9 and to a and j3. Then, the new distribution function which we denote by Fi{a, (3) 
will look like 

Fi(a, (3) = Fo(T(a, /3), P(a, P))J{a, /3), (14) 

where J{a,P) is the Jacobian 



Jia,P) 



dT dT 

da 80 

dP dP 

da dp 



(15) 



Finally, within our model, the neutrino flux observed on the Earth ($i) and the flux predicted 
by the standard solar model ($o) are related simply by 



12 



$1 ^ Fi(7r/2,7r/2) 

$0 Fo(7r/2,7r/2)- ^ ^ 

Now, we proceed to the calculation of the angular distribution function Fi. To do that, 

we need to find the defiection angles, that is the functions T{a,l3) and P(a,/3). For this 

purpose, let us consider the equations of motion for a charged neutrino in the magnetic field 

dp ee. ^^. dE ^ . ^, 

where p and v are the neutrino's momentum and velocity, respectively, ee is the neutrino 
electric charge. Since the electron neutrino mass is experimentally known to be less than 



7^9 eV [Q, to a very high accuracy the neutrino velocity equals the velocity of light, 
V ^ c. As for the neutrino energy, E, its magnitude depends on the reaction in which the 
neutrino is born. On the other hand, it must be greater than the experimental threshhold 
which is different for different experiments, ranging from about 0.2 MeV to 7 MeV. 

To simplify the equation of motion, recall that we are only interested in what happens 
in a close vicinity of the direction Sun-Earth, that is near the direction 6 = 7r/2, (p = 7r/2. 
Therefore we can neglect the curvature of the magnetic lines of force and of the convective 
zone containing these lines of force and approximate the convective zone by a flat slab of the 
same thickness, D ^ 2 x lO^^cm, containing the straight lines of force. This slab is placed 
perpendicular to the axis Earth-Sun. The distance between the centre of the Sun and the 
middle of the slab is taken to be rf ^ 6 x 10^'^cm - which is the distance between the centre of 
the Sun and the middle of the convective zone. Furthermore, we assume that the magnetic 
fleld is zero outside this slab whereas inside the slab the only nonzero component is H^ (i.e., 
we ignore the poloidal magnetic fleld altogether). In the absence of a detailed model for z- 
dependence of the magnetic fleld, we adopt the simplest form, i.e. a linear variation with z: 

Hy = H, = 0, H, = {^)z. (18) 

where (^^) is the average vertical gradient of the large-scale magnetic fleld. 

Since the deflection angles are assumed to be small, we can neglect the variation of 
magnetic fleld over the trajectory of any single neutrino flying at some 9. In other words, we 
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assume that all the neutrinos flying at the same angle 6 will feel the same strength of the 
magnetic fleld, while the neutrinos flying at different angles 6 will feel different strengths 
of the magnetic fleld. To paraphrase it once more, we assume that a neutrino flying at an 
angle 6 upon entering the convective zone will move in the uniform magnetic fleld all the 
way through until it leaves the Sun; the strength of that uniform magnetic fleld is taken to 
be the same as the strength of the magnetic fleld at the "entrance point": 

Hy = H, = 0, H, = {^)dcot 9 ^ {^)d{\ - 6). (19) 

Now, the motion of a charged particle in a uniform magnetic fleld is well-known: along 
the magnetic fleld the motion is free, while in the perpendicular plane the trajectories, both 
in coordinate and momentum spaces, are circles orbited with Larmour frequency. Thus, if 
a neutrino enters the convection zone with the speed components foxj'^oyj'^oz; it will leave 
the Sun with the speed components 

Vx = Vox, (20) 

Vy = VQy COS OUT + VqzVQz sin UJT ^ VQy + f Oj;(5, (21) 

Vz = —Voy sin LOT + Voz COS UJT ^ —Voy5 + Vqz- (22) 

Here, r = D /c is the time of flight through the convection zoneQ, 5 = lot is the angle by 
which the magnetic fleld rotates the neutrino speed vector projected to yz- plane; uj is the 
Larmour frequency: 

,7r , eec.dHx, , ._„. 

uj = ujocota^ {--a), uJo = —{——}d. (23) 



^Note that the actual times of flight, besides being diflFerent for different neutrinos, are in fact 
longer than D/c because the neutrino paths are not straight, but this difference is negligible (and 
even if it were not, assuming r = D/c for all neutrinos would only decrease the effect rather than 
increase it). 
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From these equations we see that the condition of vahdity of our small-angle approxima- 
tion can be formulated as sin 5 ^ 6. 

Having obtained the solution of the equations of motion, we are now able to find the 
relation between the initial angles 6, cj) at which a neutrino enters the convective zone and 
the final angles a, j3 at which the neutrino leaves the Sun. For this purpose, let us express 
the initial speed components v^x^v^y^voz through the initial angles 6*, 0, and the final speed 
components Vx,Vy,Vz through the final angles a,/3: 

vqx = csinO coscf), voy = c sin 9 sin cj), vqz = ccos9. (24) 

Vx = c sin a cos /9, Vy = c sin a sin /3, v^ = ccosa. (25) 

Assuming as usual that the deflection angles are small, we can write 

e = a + Aa,(j) = (3 + A(3, (26) 



where Aa, A(3 are small. Now, solving the system of equations, Eqs. ( pO] ) - (^) with respect 
to Aa, A/3, we find 

TT 

9 = a ~ 5osin/3cot « = T(«,/3) ^ a — So{— — a), (27) 

(f) = 13- 6ocot^ a cos (3 = P{a, (3)^ /3, (28) 

where 6o = uqt. Having expressed initial angles in terms of final angles, that is, having 
obtained the functions T(a,/5), P{a,l3) we can now calculate the Jacobian in Eq. ([T5|) and 
then insert the result, together with Eqs. (pTl) , (^) into our basic formula, Eq. (|^). Lastly, 
we recall Eq. ([T6|) giving the neutrino flux observed on the Earth within our model. 
Thus, finally, the neutrino flux observed on the Earth within our model is given by 

^ = 1 + ^0 = 1 + ^^. (29) 

An important thing about this result is that it is the field gradient rather than magnetic 
field itself that causes the effect. Indeed, using our method, one can show that if the toroidal 
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magnetic field in the convective zone were uniform, it would not lead to any substantial 
anisotropy of the neutrino flux at all. 

To compare Eq. ( pQ]) with the experimental data, we need to know the value of the 
mean gradient (along the solar rotation axis) of the large-scale toroidal magnetic field in 
the convective zone , {^§^)- Unfortunately, the magnitude of this gradient is not very well 
known, so let us try to reverse our problem and ask: what value of the gradient will be needed 
for our mechanism to explain the solar neutrino deficit? Assuming the neutrino charge to 
be e = 10~^'^, its energy E = 0.8 MeV and requiring the flux deficit to be 60 = —0.5, from 
Eq. ( p9D we find that the gradient must be 

(±^\ ~ _i.i X 10-^ G/cm. (30) 

oz 

Is it a reasonable figure or not? A crude estimate of the gradient can be obtained by dividing 
H hj h where H is the maximum value of the magnetic field reached at the latitudes of 
about ±10° [Q and h is the distance from that latitude to the solar equatorial plane, 
/i = t/sinlO°^ lO^^cm. 

As for the possible value oi H, it is a subject of a debated controversy. On the one 
hand, it is claimed [^ that values of H greater than lO^G are ruled out by the non-linear 



growth-limiting effects; on the other hand, there are arguments based on the helioseismology 
data that it can reach as large values as a few million G 1^^. Anyway, magnetic fields up 
to 10'* G (or even 10^ G I^S]) are widely used by many authors trying to explain the solar 



neutrino puzzle. So, we leave it to the reader to make his/her own judgement on that point. 
Note also, that it is the magnetic field close to the surface of the Sun which reaches its 
maximum at 10° latitude, and this latitude may be higher (or lower) for magnetic fields 
located at larger depths. That brings in an additional uncertainty to the estimate of the 
gradient. If we do admit that the magnetic field in the convective zone may vary in the 
range H = 10^ -^ lO^G than the value of the gradient may vary in the range 

oz h 
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Hence we see that the value of gradient needed to explain the neutrino deficit, Eq. (|30|) may 
indeed exist in the convective zone of the Sun. 

A few remarks are now in order. 

l.We have just considered the case of a = 7r/2 corresponding to the case when the Earth 
traverses the plane of the solar equator. But that happens only twice a year: in June and 
December - the days of summer and winter solstice. During the rest of time the direction 
Sun-Earth is at an angle to the solar equatorial plane. Fortunately, this angle is rather small: 
\tt/2 — «| ^ 7°, so that our small- angle approximation still applies. (The boundaries are 
reached when the Sun - Earth axis lies in the northern solar hemisphere at an angle of 7° in 
September, autumn equinox, and through the southern solar hemisphere, at the same angle, 
in March, spring equinox.) If we keep a- dependence in our calculation then we obtain that 
the deficit of the neutrino fiux varies according to 

|i = l + 4-. (32) 

$0 sm^a 

Thus for extreme values of a we have^: 

|i = 1 + -^^- = 1 + —^ ^ 1 + 1.01 X 6o. (33) 

$0 sin2 97° sin^SS" ^ ^ 

Recall that this result is valid within our system of approximations: we neglected the cur- 
vature of the solar convective zone and substituted the true value of the magnetic field 
gradient by its average value. But even if these assumptions are dropped, there are hardly 
any reasons to expect that this result would change dramatically. 

So we can conclude that the seasonal variations of the neutrino fiux predicted by the 
present model are rather small which makes them very hard to observe in the experiment. 

2. All the above derivation was of course completely classical. Indeed, one can reasonably 
expect that quantum corrections (e.g., due to neutrino diffraction) must be vanishingly 
small. 
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Note that Eq. (32) cannot be extrapolated to large angles so it does not mean that the neutrino 



flux is zero at a = vr and a = 0. 
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3. One might wonder about possible influence of other magnetic fields encountered by 
the neutrinos on their long way from the Sun to the Earth. Specifically, there are the solar 
magnetic field outside the Sun, the interplanetary magnetic field and, finally, the terrestrial 
magnetic field. 

First of all, the terrestrial magnetic field is completely transparent for solar neutrinos 
because the penetration ability of a particle of charge q and momentum p is controlled by 
the factor pc/q which is huge for a neutrino with the charge 10~^^ - it is by many orders of 
magnitude greater than this factor for the penetrating particles of cosmic rays. The neutrino 
deflection angle due to Earth's magnetism is also negligibly small since both the terrestrial 
magnetic field and the time of neutrino flight through the terrestrial magnetosphere are 
much less than those for the convective zone of the Sun. Similarly, the influence of the outer 
magnetic field of the Sun and the interplanetary magnetic field can be neglected: although 
the time of flight outside the Sun is about 10^ longer than inside, the magnetic field outside 
is 10~^ -7- 10^^ G i.e., many orders of magnitude weaker than inside, so that the product of 
these two factors outside the Sun is much smaller than inside. In addition to that, let us 
stress once again that we need large field gradients rather than large magnetic fields which 
are unlikely to emerge in the interplanetary space. 

4. Also, the role of various electric fields, both inside and outside the Sun, seems to be 
negligible for the present problem. Inside the Sun, as is well known for any plasma, any 
electric field dissipates very quickly due to high conductivity of the plasma]^. Outside the 
Sun, the electric field could exist, for instance, due to nonzero total electric charge of the 
Sun. Although there is not much information available about the strength of such electric 
field, we can be quite confident that it cannot be strong enough to capture the neutrinos 



^'^In fact, there may exist exceptions to this rule. Under certain circumstances (for example, in the 
conditions of solar flares) the electrical resistivity may be greatly increased so that a local electric 
field may arise [p5l. This point is probably worth further studying. 



with the charge 10~^^e and the energy i? ~ 1 MeV before they reach the Earth (or even 
decrease their energy to any extent). If that was possible, then the solar wind could never 
reach close to the Earth, too. The terrestrial electric field, about 1 V/cm in strength, is also 
helpless to stop the neutrinos. 

5. One might wonder if the nonzero neutrino charge would affect the solar neutrino 
detection at Kamiokande. The point is that there arises an electromagnetic contribution to 
the cross section of the neutrino-electron scattering which might modify the Kamiokande 
results. However, HQ,, = lO^^'^e then the cross-section of the electromagnetic UgC scattering 
is about 10 orders of magnitude smaller than the cross-section of weak Uge scattering and 
therefore can be completely neglected. 

6. At this stage, we have not considered the effect of solar matter on our results. 

7. Now we come to the discussion of the most serious flaw of the suggested scenario in 
its present form. The point is that up to now we tacitly assumed that the direction of the 
magnetic field is such that the value of the deficit, (— ^o); is positive. But recall that the 
large-scale magnetic field in the convective zone reverses every 11 years which means that 
the gradient, too, changes its sign every 11 years. That means that taken as it is, Eq. ( PU[ ) 
would predict that each 11 year period of neutrino flux deficiency must be followed by 11 
year period of neutrino flux excess of the same magnitude so that the flux averaged over the 
22 year cycle would be the same as predicted by the Standard Solar Model. 

One can think of several possible ways out of that difficulty. 

The most natural one is to recall that our previous calculation was based on the ap- 
proximation of small deflection angles, or, more precisely, the smallness of the parameter 
^0- We can expect that this small-angle approximation is valid as long as sm6o ~ So or 
|<^o| ~ 0.7. But if the gradient of the magnetic field is greater than 2 x 10^^ G/cm than our 
approximation does not work anymore and a more exact calculation is needed. 

Naively, one might expect that the neutrino deficiency must alternate with the neutrino 
excess at 11 year intervals independently of the magnitude of the gradient: just note that 
when the magnetic field configuration is defocusing, one would expect the neutrino deficiency 
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and when it is focusing, the neutrino excess. Each reversal of the magnetic field means a 
switch between focusing and defocusing modes so that any 11 year "deficiency" cycle would 
be followed by the 11 year "excess" cycle, however great the gradient of the magnetic field is. 
Nevertheless, there are arguments based on simple geometrical optics considerations which 
show that it is not the case and if the gradient is large enough then the neutrino deficiency 
can occur both for the defocusing and the focusing configuration! 

Another option is to try to relax the solar upper bound on the possible electric charge 
of the electron neutrino obtained in |^, since the neutrino charge and the magnetic field 



gradient come always as the product e x (^^) rather than separately. Note that if we used 
only the most reliable bound (i.e., direct experimental bound extracted from the elastic z/gC 
scattering data) on the electron neutrino charge, e ^ 3 x 10^^*^ [^, the required value of 



the gradient would be less than the above estimate by the factor of 300. 

Let us also mention briefly that at present we cannot rule out the possible existence of 
a primordial magnetic field of as much as 10^ G inside the core of the Sun [^]. Within the 
present context, it would very interesting if any evidence could be obtained concerning the 
existence of significant gradients of that field near the plane of the solar equator. 

Finally, although it is not as much appealing, we should not discard the possibility that 
our mechanism is effective only during alternative 11 cycles or even only during the periods 
of active sun within the alternative 11 year cycles while some other mechanism is responsible 
for neutrino depletion during the rest of the time. This possibility will have to be considered 
much more seriously if the anticorrelation of the neutrino deficiency with solar activity is 
established firmly by the future experiments. 

V. DISCUSSION 

As we stressed in the preceding sections, our present knowledge of the structure of the 
solar magnetic field is rather limited. Thus, we cannot rule out such values of the magnetic 
field that would lead to large deflections of the neutrinos travelling through the convective 
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zone of the sun. The quantitative theory for this case seems more difficult to construct and 
we do not attempt it here. Yet, it is instructive to discuss here some quahtative features 
of such theory based on simple physical considerations, keeping in mind the results of four 
different solar neutrino experiments available by now. 

Before doing that, let us very briefly summarize the solar neutrino data. 

1) Anticorrelation of the neutrino flux with solar activity is probably observed in the 
Homestake data [^,0 . 



2)No such anticorrelation is observed in the Kamiokande data f^. 

3)The higher neutrino flux (i.e., less neutrino deficit) is observed in Kamiokande experi- 
ment than in Homestake experiment. 

4) The higher neutrino fiux is observed in SAGE [0] and GALLEX |^ experiments than 
in Homestake experiment. 

It is very important for us to note that the experimental thresholds of neutrino energy 
are rather different in those experiments: Efjome = 0.816 MeV, Exam ~ 7.5 MeV, and 

^Ga/Zex = 0.233 McV. 

Now, let us discuss qualitatively some effects that are consequences of our hypothesis. 
These effects are controlled by the neutrino electric charge e, gradient of the solar magnetic 
field (^^), and the neutrino energy E (for the moment, we forget about other relevant 
parameters). However, these quantities enter not separately but only through the ratio 
e{^^) / E. Therefore, changing the magnetic field will be equivalent to changing the neutrino 
energy correspondingly. Furthermore, it is natural to assume that both neutrino fiux deficit 
and anticorrelations grow with the increase of the gradient. Hence we obtain that the 
anticorrelations have to be smaller for more energetic neutrinos. And this is exactly what 
is needed to qualitatively explain the difference between Homestake and Kamiokande data 
(see 1) and 2) above). Also, by the same reasoning, within our scenario one can expect less 
deficit in Kamiokande than in the Homestake experiment. Thus we can summarize that it 
is rather plausible that the present scenario can, in principle, account for the three out of 
four main experimental features, see l)-3) above. 
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Now, as for the fourth feature, i.e., results of gaUium experiments, our hypothesis seems 
to predict greater deficit than Homestake and thus looks disfavored by gallium results. How- 
ever, one must remember that: 1) the difference between Gallium and Homestake results, 
from the viewpoint of our hypothesis, must be less pronounced than the difference between 
Homestake and Kamiokande data. This follows from the fact that the ratio of the charac- 
teristic neutrino momenta for Homestake-Gallium data are smaller than for Kamiokande- 
Homestake data; 2) The errors of Gallium data are still larger than those of Homestake 
data. 

Now, we would like to draw attention to a curious coincidence in the solar neutrino 
data. Kamiokande does not see ant icorrelat ions during the whole period of its operation, 
i.e., 1987-1993 (part of solar cycle 7^ 22). And, according to |^^ there are no anticorrelation 
in Homestake data during the years 1970-1977 (part of solar cycle # 20). Also, the latest 
data do not confirm the anticorrelation: large number of the sunspots in 1991-1992 was 



accompanied by high counting rate [^. Therefore, one is tempted to speculate that, due to 
some reason, the anticorrelations are much more prominent in the odd-numbered solar cycles 
while being suppressed in the even-numhered cycles. If we take this conjecture seriously, 
it would be easy to conclude that the neutrino-depleting mechanism must somehow be 
correlated not only with the strength of the solar magnetic field but also with the direction 
of the toroidal solar magnetic field which reverses every 11 years. Obviously, this feature 
would be difficult to accomodate within any of the existing scenarios except the present one. 

VI. A SECOND NEUTRINO FLUX 

Apart from the reduction of the conventional (i.e., thermonuclear) neutrino flux, a spec- 
tacular feature of our scenario is the prediction of a "second flux" of electron neutrinos and 
antineutrinos from the Sun. While thermonuclear neutrinos are produced due to the weak 
interactions of the neutrino, the second flux arises due to the electromagnetic production 
of neutrino- antineutrino pairs. The most important process would be that of plasmon de- 
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cay into a neutrino-antineutrino pair. Thus the second flux would consist of low-energy 
(about 200 eV) neutrinos produced in plasmon decays in the core of the Sun, the number 
of such neutrinos being much greater than that of the thermonuclear neutrinos. It would be 
very interesting to consider the posibility of detecting this second neutrino flux, about 10^^ 
s'^^cm^'^ in magnitude, on the Earth. 

VII. CONCLUSION 

To conclude, in the context of the solar neutrino problem we studied the consequences 
of the hypothesis that the electron neutrino has a small but non-vanishing electric charge. 
The main general consequence is that the solar neutrino flux can be anisotropic. The cause 
of that anisotropy is the antisymmetry of the large scale toroidal solar magnetic fleld about 
the solar equatorial plane which leads to a large gradient of the magnetic field along the 
direction normal to the solar equatorial plane, denoted by (^j^)- It is this gradient that 
results in anisotropic defiection of the neutrino fiux. The magnitude of the anisotropy is 
controlled by the product of the gradient, {^^) and the neutrino charge, ee. Arguments 
based on the energetics of the Sun show that the neutrino charge must be less than e ^ 10^^^ 
2^ . Unfortunately, the value of the gradient (^^) is not very well known and, according 



to a rough estimate, may vary from 10^^ G/cm to 10~^ G/cm (or even perhaps up to 10~^ 
G/cm). We calculated, in the linear approximation, the deficit of the solar neutrino fiux 
observed on the Earth caused by the anisotropy of the neutrino fiux at the surface of the 
Sun. Assuming that the neutrino charge is equal to 10~^'^e and its energy E = 0.8 MeV we 
found the value of the gradient which is needed to obtain a 50% deficit by our mechanism: 
about 10~^ G/cm. If the neutrino charge is much less than 10~^'^e, it is unlikely to produce 
any observable effect under the action of the magnetic field of the convective zone. (If one 
considers the magnetic field of the core and radiative zone, it is not ruled out that much 
less charges are still interesting, but any quantitative conclusion on that point is difficult to 
reach in view of our poor knowledge of those magnetic fields). 
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We then discussed some attractive experimental implications of this scenario as well as 
the problems which have to be solved so that this scenario could be considered as a full- 
fledged solution to the solar neutrino puzzle. 

We would like to stress that our mechanism has two distinctive features which make it 
different from the other proposed solutions (such as MSW, "just so", spin-flip, spin- flavour 
precession and others): 

1) Energy dependence of the effect: the less is the neutrino energy, the more is the 
suppression of the neutrino flux; 

2) 11-year variation: the neutrino flux is anticorrelated with solar activity.|^ 

(In connection with these two features, we recall that: spin-flip mechanism is energy 
independent; MSW solution does not exhibit any 11-year variation.) 

The second feature of our scenario is shared by all the mechanisms employing the coupling 
of the neutrino magnetic moment to the solar magnetic field in the convective zone (involving 
or not involving the change of flavour). However, our scenario is different from all those 
mechanisms in an important way: the seasonal variations predicted by those mechanisms 
are not predicted by our model. The reason is that our effect is controlled by the gradient 
of the magnetic fleld while all the others depend on the magnetic fleld itself. 

On the other hand, our model differs very clearly from the MSW scenario or vacuum 
oscillations ("just so"). In these oscillatory scenarios the missing electron neutrinos are 
changed into neutrinos of other flavours (z/^ or z/^). In our model, no u^^ or z/,- arise to replace 
the missing z/g. This difference will be crucial for the forthcoming neutrino experiments 
involving neutral current reactions (e.g., SNO or BOREXINO). If an oscillatory scenario 
(MSW or "just so") is correct, the number of observed neutral current events will be close 
to the Standard Solar Model prediction (due to the additional contribution of z/^ and z/^). 



^^Note that such anticorrelation may be absent if the deflection of the neutrino flux is caused by 
the magnetic field of the solar core rather then that of the convective zone. 
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On the contrary, if our mechanism works, the observed number of such events will be less 
than the theoretical prediction. 

But, independently of whether this scenario survives or not in its present form, our 
arguments show that a more general problem of the possible anisotropy of the neutrino flux 
due to the interactions of the neutrino with the solar matter and electromagnetic fields is 
certainly worth further pursuing. Alternatively, the detailed observations of solar neutrinos 
in future experiments can provide a unique way of obtaining bounds on (or evidence of) the 
non-zero neutrino charge. 
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